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METHODS  OF  TABULATING  BALLISTIC  FUNCTIONS  BASED  ON 
THE  SQUARE  LAW  OF  DRAG 

ABSTRACT 

The  solution  of  the  equations  of  motion  of  a  projectile  subject  to 
the  Newtonian  Law  of  Resistance  is  expressible  in  terms  of  quadratures. 
By  applying  oertain  transformations  these  quadratures  are  now  greatly 
simplified,  and  the  functions  leading  to  the  desired  information  per¬ 
taining  to  the  trajectory  are  quickly  and  accurately  calculated.  There 
is,  then,  given  a  method  of  determination  of  the  basio  parameter  k  which 
leads  to  very  good  results,  as  shown  by  numerous  comparative  examples. 
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LIST  OF  SYMBOLS  AND  FORMULAS 


C  =  ra/id^ 

G  =  v/jJSj, 

H 

*D 

T  =  t(/3,  b 
x  =  x(/3,  b 
y  =y(0,  b 

b 

0 

g 

k 


) 

) 

) 


Ballistic  coefficient 
Resistance  function 
Relative  air  density 
Drag  coefficient 
Time  function 
Range  function 
Altitude  function 


£  (9)  +  g/(2kx^),  quantity  associated  with  the  square  law 


o  *  o 

c(£)  =  cos  e(£ ) 


Gravitational  aooeleration 
/•V6'  quantity  assooiated  with  the  square  law 


s 


s(£  )  -  sin  e(£) 


t 


Time 


v 


Speed  of  projectile 


Speed  of  projectile  at  point  of  departure,  summit,  point 
of  fall,  respectively 


x 


y 


Abscissa  (horizontal)  in  system  0(x,  y)  of  reference  de¬ 
scribing  the  motion  of  the  projectile 

Ordinate  (vertical)  in  system  0(x,  y) 


£  -  ( b  -  {  )A 
r  =  0  V2 

n  = 

9  Inolination  of  trajectory 

0(£)  Inverse  function  of£  (C) 

£  vs  £  (©)  seo^0  d  0 

P  Air  density  at  altitude  y 

"  Angle  of  fall 
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Intentionally  left  blank. 


INTRODUCTION 


The  application  of  the  Newtonian  Law  of  Resistance  to  ballistic 
investigations  has  been  found  to  give  exact  results  only  in  rather 
speoial  cases,  but  does  lead  to  good  approximations  in  many  more  oases 
of  praotioal  interest.  There  are  certain  projeotiles  whose  velocities 
are  limited  to  relatively  small  ranges,  and  whose  resistance  it  is 
difficult  to  determine  accurately.  As  pointed  out  by  L.  S.  Dederick 
and  F.  V.  Reno^,  S.  J.  Zaroodny  and  D.  N.  Brooks^,  and  others,  it  is 
for  projeotiles  such  as  the Be  that  the  square  law  of  drag  may  pro¬ 
fitably  be  employed  -  provided  one  assumes  the  product  of  air  density 
and  drag  ooeffioient  to  be  oanstont  over  the  whole  trajeotory. 

While  extensive  tables  of  exterior  ballistios  based  on  the  square 
law  of  drag  -  suoh  as,  for  example,  the  tables  of  Otto  and  Lardillon  - 
have  existed  far  a  long  time,  it  was,  nevertheless,  felt  that  an 
ao cur ate  tabulation  of  the  square  law  ballistic  functions  themselves 
would  not  only  be  useful,  but  would  also  fill  a  gap  hitherto  existing 
in  this  field.  It  is  for  these  reasons  that  the  job  was  begun,  early 
in  1945,  and  carried  out  on  the  Aberdeen  IBM  Relay  Calculators,  under 
the  direction  of  I.  J.  Sohoenberg. 

There  are  thus  now  available  six  tables,  whioh,  it  is  hoped  will 
greatly  facilitate  the  use  of  the  square  law  in  its  region  of  applica¬ 
bility*.  The  iYinctions  tabulated  are 

r «  , 

1.  £  (e)  =  J  sec0  ede. 

o 

This  table  is  the  same  as  the  one  published  in  Cranz4,  but  improved  by 
the  weeding  out  of  numerous  errors.  In  this  table,  whioh  gives  £  (9)  for 

O  '  ,  ■  a  O  1 

0  =  0  0  (1  }  87  0  ,  the  number  of  decimal  places  e^diibited  drops  from 
seven  for  low  9  to  four  for  large  values  of  9. 

2.  s(4)  =  sin  9(£),  to  nine  deoimals,  for  £  =0.00  (0f0l)  50*00. 
Here  9(0  is  the  inverse  IVnctian  of  £ (9). 

3.  c(£)  =  cos  9(£),  to  eight  decimals,  for  ihhe  same  range  of  £  * 

4.  X(/3t b)  =/^c(£)  d^/tf  . 

1 

^  Dederick  -  Reno,  “Solution  of  Differential  Equations  .of  Motion  of  a 
Projectile  with  Newtonian  Resistance  in  an  Atmosphere  srf  Variable 
Density",  BRL  No.  243,  1942. 

^  Zaroodny  -  Brooks,  "Charts  for  the  Exterior  Ballistics  of  Mortar  Fire 
Based  on  the  Square. Law  of  Drag",  BRL  No.  661,  1948, 

3 

The  Ballistic  Research  Laboratories  have  a  limited  supply  of  thede 
tables  for  distribution,  upon  request,  to  those  ageneiaa  requiring 
them. 

*  C.  Cranz,  "Lehrbuch  der  Ballistik",  Vol.  1,  Berlin,  1S25* 
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.  Y(/3 ;b)  =A)  d^. 

6*  T(/3,b)  =/*£(£)  d/0/^f? 
1 


These  last  three  functions  X,  Y,  and  T,  have  been  tabulated  to 
eight  deoimal  places,  for  b  =0,1  (0,l)  2.0,  ^=0,00  (0,02)^3,00,  but 
with  the  upper  limits  of  these  ranges  restricted  by^(b-l)  *  2,  as 
shown  in  Fig.  2,  p.  16,  To  facilitate  interpolations  all  tables  also 
show  the  differences  up  to  the  third!  order. 

Most  computations  were  conducted  to  eight  or  more  decimal  places, 
and  oheoked  by  differencing,  in  the  case  of  X,  Y,  T,  with  respect  to  b 
as  well  as /3  .  A  detailed  description  of  the  oonstruotian  of  each  table 
will  be  found  in  the  later  seotians.  There  is,  finally,  suggested  a 
method  of  obtaining  the  basio  quantity  k  fr cm  the  ballistio  ooeffioient 
C  whioh  leads  to  very  good  results  for  the  range  and  time  of  flight, 
as  shown  by  various  examples. 


THE  DIFFERENTIAL  EQUATIONS  OF  THE  TRAJECTORY 

The  motion  of  the  center  of  mass  of  a  projectile  will  in  the  fol¬ 
lowing  be  referred  to  a  system  of  Cartesian  coordinates  0  (x,y)  having 
its  origin  0  at  the  point  of  departure,  the  x-axis  being  horizontal 
and  positive  in  the  direction  of  fire,  and  the  y-axis  being  vertical 
and  positive  upwards j  see  Fig.  1. 


Fig.  1 
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In  Buch  a  system  of  coordinates  the  equations  of  motion  may  be 
written  in  the  •well-known  form 

X  e  -  Ex 

(1) 

y  -  -  EJr  -  g, 

where  the  quantities  i,  x,  eto.  ocourring  in  these  equations  have  the 
following  meanings i 

x  —  dx/dt,  t  being  the  time, 

it  =  d2x/dt2,  y  *=  dy/dt, 
g  is  the  gravitational  aooeleration,  and 

E  =  G  H/fe.  (2) 

Here 

H  =  ®"ay 

is  the  ratio  of  air  density  at  altitude  y  to  that  at  sea  level, 

8  =  Vo*!, 

is  the  resistance  funotian,  the  drag  ooeffioient,  v  the  speed  of  the 
projeotile,  and  C  the  ballistic  ooeffioient.  Making  use  of  these  formulas 
for  G,  H,  and  C,E  beoomes 
E  =  kv, 

where 

k  (3) 

Since  K^j  is  dimensionless,  k  itself  has  the  dimension  of  a  reciprocal 
length.  "While  k,  by  (3),  is  actually  dependent  on  altitude  y  and  Maoh 
number  M,  it  will  here  be  assumed  oonstant  along  the  whole  trajeotory. 

If,  under  this  assumption,  the  angle  0  of  inclination  of  the  tra¬ 
jectory  -  rather  than  the  time  t  -  is  taken  as  the  independent  variable, 
the  equations  of  motion  may  be  replaced  by  the  equivalent  system^ 

gd(v  cos  0)  =  kv3  d  0 

gdx  =  -v  d  8 

2  <4> 
gdy  —  -v  tan  0  d  0 

gdt  =  -v  seo  0  d  0. 

^  This  is  shown,  e.g.,  in  Cranz,  "Lehrbuch  der  Ballistik,  p.  110. 
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The  "principal  equation"  (4,1)*  integrates  into 

g(v  oos  ©)  2  =  -2k  £(©)  +  2kb,  (6) 

■with 

-©  * 

£(©)  =  J  seo0  ©  d  ©,  (6) 

o 

and  2kb  being  the  oonstant  of  integration.  For  the  initial  conditions l 
at  t  =  tQ,  x  =  xo,  y  =  v  =  vQ,  ©  -  ©q,  b  obviously  be o ernes 

b  =  g/(2kx2)  +£o  (7) 

with 

« O  -■<  <«0>. 

X  "  T  oos  ©  . 

O  O  0 

Sinoe  at  the  summit  S  of  the  trajeotory  ©  =  0,  v  *=  v  ,  b  may  also  be 
oharaoterized  by 

b  =  g/(2kv2).  (8) 

Solving  eq.  (5)  for  v  leads  to 
v2  *  g  sec2  ©/2k  (b  -£(©)), 
whenoe 

2  k  dx  =  -  seo2  ©  d  ©/(b  -  £  ) 

2  k  dy  =  -  ton  ©  seo2  ©  d  ©/{b  -£)  N  (®) 

(2kg)^2dt  =  -  seo2  ©  d  ©/(b 

TRANSFORMATION  OF  THE  DIFFERENTIAL  EQUATIONS 

While  the  solution  of  the  problem  is  thus  reduced  to  quadratures, 
seme  further  simplications  may  be  aohieved  by  introducing  the  inverse 
funotion  ©  (£)  of£(©).  Let  us  put 

o(£)  =  cos  ©  (O.  s(C  =  sin  ©(4)  (10) 

By  virtue  of  d  8  =o!(Od$  above  equations  (9)  take  on  the  form 

2  k  dx  »  -  eC£  )  d£/(b  -£) 

2  k  dy  =-  s(()  d£/(b-0 

(2ke)1/2dt  ^  -  o«  )d$/(t  -j  )1/2 

1  By  equation  (4.l)  is  meant  the  first  equation  of  set  (4). 

8 


(11) 


If  we  now  define  a  new  variable  by  means  of 

b  -£  =  b/3, 

i.e.  put 

£=b(l-/3),  (12) 

equation* (ll)  further  simplify  to 

2fc  dx  =  o  ( £  )  d /0//0 

2  k  dy  =  s(£  )  */3//3  (13) 

2  lev  dt  -  c(£  )  d/3//(3 

By  virtue  of  equations  (5)  and  (8)  the  variable  turnB  out  to  be 

^=r  Vg  aeo^  «/v2«  (14) 

Therefore,  as  deoreases  from  Q  to  $  *=  0,  <£  decreases  from  £  Q  to  0, 

and /G  increases  from 

/30  =  (vs  seo  Qjrf  =  1  -  £  Jb 

to  /2>  —  1;  as  the  projectile  moves  down  the  descending  branch /Q  is  in¬ 

creased  further.  In  terms  of  the  parameter the  trajeotory  is  now 
represented  by 

2k(x  -  xQ)  =X(/3,b)  -  X(/3q»  b) 

2k (y  -  y0)  =  n/3, b)  -  Y(/90,b)  (15) 

2kvB(t  -  to)  =T(^(b)  -  T (/30,  b) 

where 

z(/S,b)  =f%(£)A/3//S 

1 

T(A  b>  (16) 

K/3,  b)  =^(4)a p/p'-b 

and 

£  ==  b(l  -/ 3 ). 

The  integrations,  desoribed  in  the  following  sections,  were  actually 
carried  out  not  on  the  integrals  (16),  but  on  related  integrals,  obtained 
as  follows i  First  it  is  advantageous  to  replace  £  by 

7}  =  -£  =b(£-  1).  (17) 
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Since  ©(£)  is  an  odd  function,  s(£)  is  also  odd,  while  o(£)  is  even* 
*(£)  =  -  b{v));  c(£)=o(ij). 

Further  we  may  write 

X(£,b)  =  X(  0f  h)  +  o(b)ln£,  (18) 

where 

*(&b)  -  o(b)J  d/8/,3,  (19) 

and 

Y(£b)  =  Y(#b)  +  s(b)  In  G  (20) 

with 

Y(#b)  =^-3(7)  -  s(b)}  d/3/0  •  (21) 

The  transformation  (17)  ohanges  the  time  integral  (16,3)  into 

T(/3,b)  =/VA 

It  is  convenient  here  to  introduce 

y  =  /3l/2. 


so  that 

T(/3,b)  =/^(7,  b)dy  (22) 

with 


S(7,  b)  =  2c(7J  ),  ?)=b(y2-l),  y2  ~  0  . 

The  problem  is  thus  reduced  to  the  calculation  of  X,  Y,  and  T# 

As  the  basis  for  all  neoessary  integrations  there  was  used  the 
formula 

A'd  j3  =  (h/2)  (f  +  f  )  +  (h2/io)  (}o  -  K)  + 

po 

(h3/L20)  (f*  +  f*)  +  R,  (23) 

where  dots  (*)  denote  differentiation  with  respect  to  jQ  .  This  formula 
may  be  derived  as  follows  1  Expanding  A  £  in  Taylor  series  about  /3 Q* 
/3  ^  +  h,  respectively,  we  have 

A*  =  ?  (hVii)  f  ^ 
o  i=l  '  0 
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and 

A° =  Si  t(-^i"lhi/id fi(1  1)# 

By  addition,  then,  ' 

A’o  =  (h/2)  (f  +  fx)  +  (h2/4)  (f0  -  fa)  +  (h3A2)  (fQ  +  fj)  +  0(h4). 

2  2 
If  we  now  split  up  the  h  term  into  two  parts  one  of  whioh  Is  (h  / 10 )  x 

(fQ  -  f-),  and  expand  the  other  part,  again  by  a  Taylor  series,  we  get 

the  formula  (23)  with 

R  =  0(h4). 

Having  obtained  the  X,  Y,  and  T  for  certain  equidistant  key  arguments 
a  subtabulation  was  performed  by  means  of  second  order  osoulation, 
in  the  following  manner!  In  order  to  calculate  the  value  of  a  function 
g(t)  tabulated  for  equidistant  t  =  t^  by  means  of  osoulatory  inter¬ 
polation  we  put 

t  *  tQ  +Th,  h  -  tx  -  tQ 

g(to+  rh)  =^(r)  (24) 

and  get* 

^(T)=D  +  D1T+..,-tDcT5  (25) 

9  oi  b 

where  the  D.  are  the  following  functions  of  _  <^(j),  j  B  0,  1, 

1  j  “  ' 

<f> j  » 

Do=*o 

B1  =  #l=hBo’ 

D2  =  (l/2)£*=  (h2/2)gy. 

®3  *  -10  4>0  -  *  1°^!  '  +  (lA)^i'  -  (3/2)^y 

D4  =  15  <j>o  +  B<t>'o  -  15<^1  +  7  <)>  ’  -  (1/2)  </>  *'  +  (3/2)^^' 

D5  =“6^o  "  1+  6^1  “  +  *  (1/2)0  o'- 

Rearranging  the  terms  in  formula  (25)  so  as  to  exhibit  the  factors  of 
q,  <£  eto.,  there  results 

4>(r)  =  <£0  p0(t)  +  <*>'  P^T)  +  </>  ”p2(t)  10.o(r)  +4>’Q1(t) 

i"  Qg  (r)  (26) 

*  For  a  derivation  of  this  formula,  see,  for  example,  LOIKIN,  M. , 

"Inversion  on  the  Eniac  Using  Osoulatory  Interpolation?  BRL  632,  1947. 
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■with 


=  1  -  lor3  +  15T4  -  6T6 

pl(r> 

*=  T(1  -  6T2  + 

8 r3  -  3r4) 

P2(r) 

=  (1/2)  r 2  (1 

-T)3 

Q0(r) 

=  l-Po(r) 

V*) 

«=  -t3(4  -  7  r  « 

*  3T2) 

Vr) 

=  (1/2)  r3  (1  - 

r)2. 

Poor  subtabulation  to  AT  e  0.1  these  polynomials  P^,  assume  the 
values  shown  in  Table  I.  * 

TABLE  X.  INTERPOLATION  COEFFICIENTS 


P0(r) 

Pi<T) 

P2(r) 

Q0(r) 

^(r) 

«2(t-) 

0.1 

0.991440 

0.094770 

0.003645 

0.008660 

-0.003330 

0.000406 

0.2 

0.942080 

0.163840 

0.010240 

0.067920 

-0.021760 

0.002560 

0.3 

0.836920 

0.195510 

0.01543  5 

0.163080 

-Q. 068690 

0.006615 

0.4 

0.682560 

0.190080 

0.017280 

0.317440 

-0.107520 

0.011620 

0.5 

0. 500000 

0.156250 

0.015625 

0. 500000 

-0.156250 

0.015625 

0.6 

0.317440 

0.107520 

0.011520 

0.682560 

-0.190080 

0.017280 

0.7 

0.163080 

0.058590 

0.006615 

0.836920 

-0.195510 

0.015435 

0.8 

0.057920 

0.021760 

0.002560 

0.942080 

-0.163840 

0.010240 

w 

0.008560 

0.003330 

0,000406 

0.991440 

-0.094770 

0.003645 

'  COMPUTATION  OF  THE  TTJNCTION  £  («) 

Since  £  (0),  as  defined  by  (6),  is  an  odd  funotian  of  ©,  it  is 
positive  along  the  ascending  branch,  vanishes  at  the  summit,  and  be¬ 
comes  negative  along  the  descending  branch  of  the  trajeotory.  It  may 
be  written  explicitly  as 

£  (®)  =  (l/2)  ftan  ®  sec  ®  +  log  nat  tan.  (-IT A  +  ®/2)^»  (27) 

£ (0)  has  been  tabulated  in  Cranz,  for^every  minute  of  arc,  from  0=0° 
to  ©  «=  87°0’,  to  8  or  less  significant  figures.  Upon  differenoing  this 
table  up  to  differences  of  third  order  two  types  of  errors  showed  up, due, 
obviously,  either  to  misprints  or  small  miscalculations.  A  listing  of 
these  values  of  £  (©)  that  were  found  to  be  in  error  by  6  or  more  in  the 
last  decimal  shown  is  given  in  Table  II.  The  errors,  when  not  due  to 
misprint,  were  eliminated  by  re computation  on  the  basis  of  eq.  (27),  in 
the  following  manner  1  Making  use  of  the  nine  plaoe  UPA  Table  of  tri- 
gancsnetrio  functions  for  each  0.0001  radians  of  argument,  we  write.. 
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TABLE  II.  ERRORS  IN  CRANZ'S  TABLE  OF  £  (©) 


Q 

CRANZ 

CORRECT 

VALCTE 

ERROR 

30 

0.056 

2500 

0.056 

2300 

+ 

200 

7 

■ 

0.131 

3  552 

0.131 

4352 

- 

800 

14 

na 

0.256 

3625 

0.256 

3525 

+ 

100 

15 

17 

0.276 

6189 

0.276 

6199 

- 

10 

17 

47 

0.326 

1715 

0.326 

1615 

100 

18 

25 

0.339 

0377 

0.339 

0337 

+ 

40 

20 

57 

0.392 

0202 

0.392 

0212 

- 

10 

21 

25 

0.402 

0083 

0.402 

0683 

- 

600 

21 

37 

0.406 

4046 

0.406 

4036 

+ 

10 

23 

55 

0.457 

6223 

0.457 

6233 

10 

25 

16 

0.488 

0669 

0.488 

9669 

9000 

25 

34 

0.496 

0831 

0.496 

0731 

+ 

100 

27 

54 

0.553 

2731 

0.553 

2631 

100 

29 

21 

0.590 

6981 

0.590 

6891 

+ 

90 

35 

42 

0.777 

3334 

0.776 

3334 

+ 

10000 

58 

27 

2.198 

352 

2.188 

352 

+1000QX 

63 

06 

2.893 

664 

2.893 

654 

+ 

1CX 

66 

39 

3.700 

597 

3.710 

597 

-1000QX 

69 

24 

4.643 

374 

4.633 

374 

41000CK 

73 

55 

7.238 

203 

7.238 

293 

- 

9  ax 

81 

03 

21.681 

00 

21.681 

08 

8XX 

83 

00 

34.811 

31 

34.811 

36 

- 

5XX 

84 

16 

51.348 

13 

51.348 

02 

11XX 

85 

23 

78.532 

3 

78.534 

1 

- 

18XXX 

86 

59' 

182.103 

0 

182.103 

7 

- 

7XXX 

0  =  0q  +  n,  ijr  =  7T/4  +  e/2  =  ifr^  +  p, 

■where  ©o,  ^  q  represent  the  first  four  deoimals  only  of  0  and  yjr, 
expressed  in  radians  to,  say,  9  decimals.  The  values  of  sin  0,  oos  0  are 
then  obtained  with  sufficient  accuracy  ft* ora 

2 

sin  0  =  (1  -  (n  /z))  sin  9q  +  n  cos  0^ 

2 

cos  ©  =  (1  -  (n  /2))  cos  0q  -  n  sin  0q 


since  the  error  ipourred  in  cutting  off  the 
shown  here  is  B  =  n^/6J ,  n  =  1.10"^. 


series  after  the  terms 


Having  thus  found  sin  0,  oos  ©,  sin  iff ,  cos  ijT ,  tan  0  seo  0  and  tan 
were  easily  oamputed.  Putting,  as  above. 


tan  )jf  =?  (tan  iff  )  +  q, 

where  (tan  lj/  )  represents  the  first  four  deoimals  only  of  tan  ^  ,  and 
using  the  notation  r  =  q/(tan  )Q,  so  that 

tan  Tjf  =  (tan  ^)0  (l  +  r)» 

In  tan  was  oaloulated  sufficiently  aoourately  from 
In  tan  \j/  =  In  (tan  )0  +  r  -  r2/2 

with  the  nine  plaoe  values  of  In  (tan  jt  ^extracted  from  the  TOA  Table, 


THE  TABLES  OF  s(£  )  and  c(^) 

If  we  consider  0  as  a  function  of  £  ,  eq.  (27)  may  be  expressed  as 
4$  =2a(£  )/c2(0  +  In  tan2  £lf/4  +  0(O/2j* 
so  that 

<£  =  Q(a)  =  (1/4)  £(1  -  s)’1  -  (l  +  s)"1  +  In  (1  +  s)  - 


From  a  table  of  Q  (s)  the  function  8(0  may  thus  be  computed 
interpolation. 

To  construct  an  eight  place  table  of  Q(s)  by  means  of  eq.  (28)  for 

the  arguments  s  »  0,0001  (O.OOOl)  0,9960  it  was  neoessary  first  to 

calculate  reciprocals  w  ~  l/u  for  u  =  0,0040  (O.OOOl)  0,9999  to  eight 

deoimals.  Using  Barlow's  Table  to  get  a  first  approximation  wq  of  w, 

correct  to  five  decimals,  a  correction  dwQ  was  obtained  in  the  usual 

manner?  If  Z  (u,  w)  =  0,  and  w  =  w  +  dw  ,  then 

o  o 

dw0  =  -  H  (u,  wo)/2w(u,  w0),  H  w  =  a  H/dw.  (29) 


In  (1  +  sj)  . 
(28) 

by  inverse 
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Since  here  Z  (u,  w)  =  u  -  w"1,  dw  =  wo(l  -  uw  ).  The  values  of  the 
natural  logarithms  needed  were  a^ain  extracted  from  the  THPA  Table. 

Having  made  this  table  an  approximate  value  s  of  s  was  then 
found  by  inversion  for  each =  0.1  (O.l)  2.0  (O.oS)  60.00  and  improved 
by 

dsQ  =  -  s0)/Za«,  sQ); 

in  this  case  Z  (£,  s)  =£  -  Q(s),  so  that 
ds0  -  -  Q(so0  /Q'(so), 

and 

Q*(s)  =  (1/4)  {(I-  -  s)'1  [l  +  (1  -  s)"1]+  (1  +  s)"1  Q  +  (1  +  s)-1]} 

This  process  was  continued  until  s  had  been  determined  aoourately  to 
nine  plaoes,  for  the  key  arguments  ^listed  above*  Next  a  subtabulation 
to  tenths  was  performed,  and  the  s{£)  obtained  f or  <£  =  0,00  (0.01)  2,00 
by  means  of  the  doubly  osculating  cubio 

£(r)  =eq  +  E1r+E2r2  +  e  3t3 

with  B  -  <t> 
o  o 

s2  =-  »*0  +  z*i  -  -ti 

E3=2^o-2^1 

giving  an  eight  place  accuraoy  in  the  interpolated  values,  with  a 
possibly  wrong  round-off  in  the  ninth  place.  The  s(£)  table  was  then 
checked  by  forming  the  advancing  differences  up  to  the  fifth  order. 

Later  calculations  required  knowledge  of  derivatives  s’  =  dx/d<£, 
s”,  and  s,n«  By  Bq.  (10) 

s'  =c4(4)  =  (1  -  S2)2. 

Therefore, 

sn  =  -  4s(l  -  s2)3, 

and 

sn«  =  4(1  -  s2)4  [6  -  7(1  -  s2) J. 

These  derivatives,  then,  could  be  computed  easily  onoe  s  =  s(0  was 
known. 

The  determination  of  the  function  o  <o.  too,  was  based  on  that  of 
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s(£  ),  s-ince 

.(«>  =£  -  .2($)J  l/2 

The  eight  place  values  of  the  square  roots  were  obtained  by  seoand 
order  interpolation 

£(r)  ^  4>0  +  (1/2)  (Ao  +  Ax)  t  +  (1/2)  A  Jr2 

in  an  eight  plaoe  table  of  roots  for  arguments  0.0001  (O.OOOl)  0.9999* 

The  derivatives  of  c,  also  needed  later,  were  obtained  by  repeated 
differentiation  of 

o?G£ )  +  s2(£  )  =  ij 

expressed  in  terras  of  the  derivatives  s'  and  s"  they  are t 

o'  -  -  ss’/° 

o"  —  -  (o'2  -f  s,2  +  ss")/0 
o'"  =  -  [3(cV  +  sf,s"  )  +  ss"  *j|  /o 
or 

o'"  =  sc7  (28o2  -15), 

where  the  last  relationship  may  be  obtained  by  differentiation  from 
,  3.4 

THE  FUNCTION  X(yg  ,  b) 

From  Eq.  (I6*l)  it  is  obvious  that  X^  o  for  /3  $  1,  for  all  values 
of  b.  As  already  pointed  out  previously,  the  actual  calculation  of  X 
was  based  on  Eqs,  ( 18)  and  (19).  The  funotion 'x(^,  b)  was  first  ocra- 
puted  for  b  =0.1  (O.l)  2.0  and  /3  =0.0  (O.l)  3.0,  with  the  upper 
ranges  of  t  restricted  by  j3  (1#-  l)  *  2.0,  as  shown  in  Figure  2. 


As  le  easily  seen  in  this  region  also  X  ^  o  for 

In  using  Eq.  (23)  for  the  integration  of  X  we  put  f(/9)  =  A(/3)  + 
B(/3)  with  M /Q)  ~  c  ( 7))/ /3  »  B (/3)  =  -  c(b )//3,  thus  considering  b 
as  a  fix;ed  parameter.  The  derivatives  of  f  occurring  in  Eq.  (23)  were 
then  obtained  from 

A  =  (be*  -  A )/jS 

A  =  (la2  c"  -  2A)//3. 

B  =  -  B//3 

3  =  -  ZB//3 


The  values  of 

»lf(o)/a^i  =diA(o.)/d/<3  1  +»iB(o)/a/5i,  i  ■  1.  2, 

which  were  needed  to  get  started,  were  found  by  L'Hospital'g  rulei 

fQ  =  -  bc'(b). 
i0  =  (b2/2)  c"0=) 

",  =  -  (b3/3>"'(b). 

By  determining  the  increment  A  on  the  basis  of  Eq.  (23),  and 

taking  advantage  of  the  fact  that 

1(1,  b)  =  0  for  all  b. 


the  function,  b)  was  oalculated  to  ten  places,  for  the  range  of 

/3  and  b  shown  above. 

The  values  of  the  increments  A  X  thus  computed  were  then  checked 

by  comparing  A  ^+2  =  A  *+1  +  A  ^+2  with  A  I*1  obtained  by  applying 
Eq,  (23)  to  hf  =  2hj  this  check  showed  the  A?,  to  be  accurate  in  most 
cages  to  eight  decimals.  The  integrations  were  checked  further  by 
differencing  the  X*s  listed  for  all  fe  ,  but  fixed  jQ  • 

These  I's  were  next  subtabulated  to  fifths  in  the  argument  , 
by  means  of  formula  (26),  which  becomes,  in  this  case 

X(/3  +  rh)  =  xqpo(t)  +  h(A0  +  V  pi(r)  +  b2  K  +  Po>  p2(t) 

+  2^  Qo(T)  +  h(A;  +  B1)  (^(r)  +  h2(A1  +  Qg(r). 

With  these  values  of  X  on  hand,  it  was,  finally,  an  easy  matter  to  get 
the  X's  themselves. 
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THE  FUNCTION  Y( ,  b) 

As  regards  the  function  Y (/3,  b)  it  is  again  apparent  from  Eq.  (16*2) 
that  Y  <  o  for  all  fc  ,  /3 .  Now  for  o  <  /3  <  1  £  «  T>  (1  -/})  <  b,  so 
that  s(£)  <  s(fc  ),  and,  therefore,  7  >■  o.  For  >  1,  on  the  other  hand, 
it  follows  from  Y_*=.7  +  s(b)  In  ,  since  Y  <  o  and  s(fc}  In  /&  >  o, 
that  necessarily  Y  <  o. 

The  region  of  integration  was  the  same  as  used  for  the  X.  In  Eq. 

(23)  we  had  to  put  here,  however, 

f(^)  =M(/3)  +  N(£) 

with 

M(/3)  =-  aiy)/,#  K(/3)  ~~  fl(b  )/& 

Consequently,. 

M  =  -  (bs*  +  M )//3  M  =  -  (fc2s"  +  ZU)//S 

N  =  -  N /£  N  =  -  2N /fi 

F  or  J3  =  0 

f Q  =  -  be* (b) 
fo  =  (b2/2)  a"(b) 

?„  =-  Cb 3A) 

The  determination  of  the  ^  7,  their  ohecking,  and  the  subsequent  sub¬ 
tabulation  were  oarried  out  in  the  same  manner  as  previously  desoribed 
for  the  function  X. 


THE  FUNCTION  T  (/?.  b) 

The  time  integral  (16.3)  behaves  similarly  as  the  range  integral 
Xf  so  that  T  ^  0  for  1.  As  mentioned  above,  it  was  oaloulated 

by  means  of  Eq.  (22).  To  make  use  of  the  integration  formula  (23)  it 
was  first  neoessary  to  determine  the  quantities 

V  -  A  V2  _  A  l/2 

hi'/l3l  &  1-1.  1=1(1)30, 

this  was  done  to  ten  decimals.  The  derivatives  of  S  with  respeot  to  Y 
are 

s«  =4  b/3  1/2  o’ (7) 

S"  =  4  b  c  *  (  77)  +  8  b  2^c "  (  7  )  . 

Checking  of  the  values  of  A  T  was  performed  in  the  same  manner  as  was 


18 


done  for  the  ^  X  and  A  Y,  and  showed  the  A  T,  too,  to  be  accurate 
generally  to  eight  decimals. 

Since  the  h  =  ^  y  depend  on  ^3  ,  the  subtabulation  of  X  to  fifths 
was  a  rather  tedious  job.  It  was,  obviously,,  necessary  to  get  for  each 
of  the  four  values  of =  C^+t2j/L00)1  ' \  0  =  0-0  (0-1)  3.0, 
j  =  1,2, 3, 4,  the  value  of  T  .  =  (  “  X  jA1*  "then  compute  the  inter- 

polational  polynomials  P.  (T^),  QAV-i)*  =  °»  1,  2,  and,  finally,  to 
use  formula  (26)  in  the  rormd  ■ 

T(/+  T^h)  =topo(tj.)  +  hs^cr^.)  +  h2S*  p2(Tj)  + 

T1  VTj>  +  hSlQl(Tj)  +  h2si  Og  ( Tj). 

TRAJECTORIES  FOR  LARGE  6  AND  C 

o 

The  tables  described  in  the  previous  sections  permit  the  calcula¬ 
tion  of  trajectories  to  a  certain  limited  extent.  For  such  angles  0  , 
or  ballistic  coefficients  C  that  cause  b,  Eq.  (7)  to  fall  outside  of0 
the  range  of  the  tables  the  trajectories  may  be  easily  obtained  by  means 
of  three  simple  quadratures,  aB  follows*  Let  us  transform  the  integrals 
(l6)  by  means  of  Eq.  (17)  to 

b>  =4  +  b) 

-y1(^,  b)  =  <*?/(? +  b)  (30) 

and 

1>  "1/2T1(^,  b)  =  "h  “1//2-/*^(^)  d?/(?  +  b)1/2. 


whence 

2k (x  -  xD)  =  X1(>7,  b)  -  X1(  ??o,  b) 

-2k (y  -  yQ)  =  \{7},  b)  -  Yx(  7 7q,  b)  (31) 

2kvs  bX/2(t  -  t^T^,  b)  -  Tx(^o,  b  )• 

The  tables  of  the  functions  c  (?)  and  s  (?)  have  been  computed  suffi/r' 
ciently  far  so  ae  to  make  the  calculation  of  the  integrals  (30)  an  easy 
matter. 

In  order  to  get  the  range  by  means  of  formulae  (3l)  one  first 
determines  ?o  and  then  computes  Y^  Ok  »*?>•.  Starting 

again  from  T}  v  O.the  quantity  ?«  >  0  is  next  calculateu  by  integrating 
to  that  value  ?  =  that  will  make  7)^,  b  )  «=  Y^(  7}  ^  b  )* 

Finally  the  range  is  obtained  by  computing  X^(?Q*'b  )  and  X^J^^  b  ), 
and  applying  formula  (31«l). 
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USE  OF  THE  TABIES 


The  proper  use  of  the  tables  described  in  the  previous  sections 
will  now  be  illustrated  by  means  of  a  few  examples;  for  the  sake  of 
simplicity  let  =  yQ  =  t  =  o. 

1.  Given  k,  0Q,  v  .  To  find  the  angle  of  fall  ,  range  x^# 
time  of  flight  t  ^  ,  an$  striking  velocity  v^ . 

p  p 

(a)  Compute  b  (0Q)  4  g/(2k  v^  cos  0Q). 

(b)  Get  4  0  =  1  -  £  qA* 

(o)  By  interpolation  find  /3  -  0^  bo  that 
0  =  2ky^  =  Y(  b)  -  Y(/?o,  b). 

(d)  By  interpolation  find  to  so  that 

-b(l 

(e)  Obtain  the  horizontal  range  xw  ft’om 

2tocw  =  X^,  b)  -  X(/3q,  b).  (31) 

(f)  Calculate  the  time  of  flight  t  w  from 

(ZkgA)1/2  tw  =  T(^,  b)  -  T(/30,  b). 

(g)  Finally,  the  striking  velocity  is  found  ft-om 

2  2 

=  g  seo  0>/2kb^w. 

2.  Given  vQ,  0  ,  x^.  This  case  can  be  reduoed  to  the  previous 
one  by  assuming  a  trial  value. of  k,  and  then  computing  the  quantities 
b,  / 3  o’  and  /&o>  *  This  is  done  until  k  satisfies  equation  (3l). 

The  quantities  tw  and  v^,  are  then  obtained  as  above. 

3.  Given  k,  vQ,  xw  .  Assume  a  trial  value  of  0O,  compute  b,  jQ  o, 

0^.  The  correct  ©Q  must  again  satisfy  equation  (3l).  L 


DETERMINATION  OF  k 

The  calculation  of  that  value  of  k  from  the  data  Cf  0  ,  v which 

will  produce  the  same  range  and  time  of  flight  as  the  one  obtained  by 

numerical  integration  is  a  difficult  job.  It  is,  however,  not  quite  so 

complicated  to  get  a  k  that  will  give  a  range  good  enough  to  be  used  as 

an  approximation.  Such  a  value  of  k  may  be  based  on  initial  and  sum- 

mital  velocity  v  and  v  ,  respectively,  as  follows i 
o  s 
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(a)  Compute  kQ  =G(vo)/VoC. 

(b)  Get  bQ  -  g/iZIcj?)  +  <£  (®0). 

1  /z 

(c)  Determine  vQs  =  (g/2k^bo)  • 

(d)  Calculate  v  =  (l/s)  (v  +  v  ). 

m  o  os 

(e)  Put  k  &  G(v^/v^p. 

2 

(f)  Compute  b  =  g/(2kxQ)  +  <£  (©0)»  and  proceed  ae  in  example  1. 

This  procedure  for  getting  k  is  obviously  justified  since,  as  seen  in 
Fig*  3,  the  motion  of  the  projectile  along  its  trajeotary  is  dupli¬ 
cated  along  the  K*  -  curve  ob  a  motion,  starting  at  0,  in  the  direction 
of  decreasing  v,  having  a  point  of  reversal  of  direotion  at  S’,  which 
is  located  slightly  past  the  summit  3,  and  terminating  at  T, 


In  Table  III  there  are  compared  summital  velocities  v  ,  ranges  x ^  and 
times  of  flight  toe  obtained  by  means  of  the  square  law  with  those  cal¬ 
culated  on  the  basis  of  the  resistance  functions  G.,  In  oases  3,4,  5 
and  6  the  values  of  v  ,  t  ,  and  x ^  ,  for  the  square  law,  were  secured 
by  means  of  formulae  ^30). 

As  the  table  shows,  the  agreement  between  the  trajectories  1b  re¬ 
markable,  the  discrepancy  amounting,  at  worst,  to  about  2%  for  the  oases 
exhibited. 
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TABLE  III,  COMPARISON  OF  TRAJECTORIES 


1 

9 

0 

(degs) 

V 

0 

(ft/sec) 

C 

v8(ft/seo) 

tw  ( 

boo) 

(yards) 

$i 

SQ.L. 

Gi 

SQ.L, 

Gi 

30 

900 

C,  =0.5 

457 

459 

22.3 

22.2 

3462 

3484 

40 

800 

=  0.6 

388 

391 

26.4 

26.4 

3475 

3502 

40 

412 

=  0.976 

281 

281 

1-5.8 

15.8 

1473 

1479 

68 

747 

c:  =  1.65 

227 

231 

39.2 

39.6 

2972 

3035 

5 

06o 

900 

C t.  «=  3.7 

340 

345 

49.2 

48.6 

54  69 

5691 

6* 

76°8V8 

639.9 

C  =  1.56 

145 

146 

36.4 

36.6 

1762 

1773 

* 

This  example  is  a  particular  oase  pertaining  to  the  105nnn  Howitaer 
far  whioh  a  epeoial  G  function  ia  available. 


'Jtt. 

Mo  Lotkin 
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